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We investigate low-lying fermion modes in SU (2) gauge theory at temperatures above the phase 
transition. Both staggered and overlap spectra reveal transitions from chaotic (random matrix) 
to integrable (Poissonian) behavior accompanied by an increasing localization of the eigenmodes. 
We show that the latter are trapped by local Polyakov loop fluctuations. Islands of such "wrong" 
Polyakov loops can therefore be viewed as defects leading to Anderson localization in gauge theories. 
We find strong similarities in the spatial profile of these localized staggered and overlap eigenmodes. 
We discuss possible interpretations of this finding and present a sparse random matrix model that 
reproduces these features. 



I. INTRODUCTION 

The vacuum of Quantum Chromodynamics (QCD) is 
a prominent nonperturbative system, whose strongly in- 
teracting nature persists even above the transition to the 
quark-gluon plasma. To get insight into its mechanism, 
spectral properties of the QCD Dirac operator are very 
useful. A nonzero density of eigenvalues at zero gives 
rise to chiral symmetry breaking via the Banks-Casher 
formula [T]. Moreover, exact zero modes are related to 
the topological charge via index theorems. 

In recent years, localization properties of the Dirac 
eigenmodes have attracted attention as they can be used 
to draw analogies to condensed matter phenomena: con- 
cepts like the mobility edge and Anderson localization 
can be studied in QCD lattice simulations. In this spirit 
the chiral transition at finite temperature^ has been con- 
jectured to be an Anderson (metal-insulator) transition. 

This goes hand in hand with different random matrix 
theory (RMT) descriptions of the Dirac spectra. In the 
low temperature phase, the existence of the chiral con- 
densate connects QCD to chiral perturbation theory and 
random matrix theory (becoming exact in the epsilon- 
regime), which explains the statistics of the low lying 
part of the Dirac spectrum [3]. 

The spectral gap in the high temperature phase, on the 
other hand, seems to call for the "soft edge" description of 
RMT, which, however, could not be supported by lattice 
data [1H7]. Instead, a transition to independent eigen- 
modes obeying Poisson statistics has been suggested [S]. 
A refined analysis by one of us has shown, that the bulk 
of the spectrum is still delocalized and subject to RMT, 
while the lowest lying eigenmodes display a transition to 
localization and Poissonian behavior of the eigenvalues 
O [To] . This effect has been confirmed to be universal 
in the sense that it does not depend on the resolution of 
the lattice. Observables, like e.g. the number of local- 



At zero temperature the situation is not clear due to the contin- 
uum limit, see |2] and references therein. 



ized modes, rather scale with the physical volume. On 
the other hand, Ref. TT] found that localization is es- 
sentially a finite volume artifact. However, they used 
a different, less restrictive definition of localization and 
thus their results are not in conflict with the rest of the 
above cited literature. 

In this work we give another crucial ingredient of An- 
derson localization in QCD, namely we identify the "de- 
fects" causing it. We show that at high temperature local 
Polyakov loop fluctuations trap low-lying modes. The av- 
erage Polyakov loop as an order parameter of the decon- 
finement (or center) phase transition approaches 1 with 
increasing temperature. Locally, however, the Polyakov 
loop takes on other values, in particular close to other 
center elements, —1 in the case of gauge group SU{2). 
The phase transition can actually be viewed as a per- 
colation of the physical sector embracing such islands 
where the Polyakov loop is close to other center elements 

The locations of these Polyakov loop fluctuations 
(which are similar to Weiss domains in ferromagnetism) 
are correlated with maxima in the profile of low-lying 
Dirac modes. We will demonstrate this for eigenmodes 
of both the overlap and staggered operator. Similarities 
between the spectrum of the overlap and staggered Dirac 
operator have already been found in the Schwinger model 
[13] and also in QCD [IS]. Our present study is, however, 
the first one when similarity of staggered and overlap 
Dirac eigenmodes is seen in lattice simulations and we 
consider this an important side result of our study. 

The observed localization can be understood via 
Polyakov loops compensating the twist caused by the 
antiperiodic boundary conditions in the temporal direc- 
tion. The corresponding Matsubara frequency is effec- 
tively lowered (locally) which results in lower eigenval- 
ues. Actually, our finding has been inspired by a similar 
localization effect in the spectrum of the gauge-covariant 
Laplace operator [TB] . The Laplacian is the square of the 
Dirac operator in the free case, so the twist picture ap- 
plies. Otherwise this operator does not share important 
chiral features like topological zero modes and conden- 
sates. 



2 



Likewise, our finding is consistent with the existence 
of a chiral condensate in the Polyakov loop sector close 
to other center elements [TTHH] (also needed for center 
symmetry breaking [25]), which implies low Dirac eigen- 
values in islands of such Polyakov loops. 

The connection of these islands to topological excita- 
tions like magnetic monopoles is attractive, but in its 
naive form contradicts the observed topological suscepti- 
bility quantitatively, see below. 

For the construction of random matrix models valid at 
high temperature we investigate the distribution of local 
Polyakov loops and find them to be uncorrelated to a 
good approximation. Hence, Polyakov loops in fact pro- 
vide the Poissonian ingredient for the Dirac spectrum. 
This is built into a novel Anderson-like random matrix 
model through supplementing it by random matrix en- 
tries that represent nearest neighbor hoppings in three- 
dimensional space. We motivate this model and show, 
that with a few parameters it reproduces the main fea- 
tures of the Dirac modes: chirality, spectral gap, RMT- 
Poisson transition and localization (to the analogue of 
local Polyakov loops). 

Our findings are based on quenched lattice simulations 
with the SU(2) gauge group, we strongly believe that 
the described phenomena are present in more realistic 
gauge theories, too. 

The paper is organized as follows. In the next sec- 
tion we describe the Dirac spectra at high temperature 
including the RMT-Poisson (chaotic to integrable) tran- 
sition and the similarity of the staggered and overlap 
modes. Sect. IIIII is devoted to the connection between 
local Polyakov loops and low-lying modes. In Sect. |IV] 
thereafter we investigate two possible interpretations of 
this finding, effective Matsubara frequencies and topolog- 
ical objects. In Section [V] we introduce and explore our 
random matrix model and finally Section [Vl| contains our 
conclusions. 



II. DIRAC SPECTRA AT HIGH 
TEMPERATURE 

We analyze quenched SU(2) lattice configurations gen- 
erated with Wilson action on a 24^ ■ 4 lattice at f3 = 2.6 
which amounts to a temperature of 2.6Tc. The aver- 
age Polyakov loop of 0.37 signals deconfinement (by 
Polyakov loop we refer to the trace of the products of 
all temporal links L{x) = 1/2 • Tr n^*=i C^o(a;o, a;) and 
we selected the physical sector of positive Polyakov loops 
by hand). 

We measured the 256 lowest eigenvalues with positive 
imaginary parts of the overlap [26, 22i (with parameter 
s = 0.4 cf. [H]) and staggered Dirac operator on 1136 
(overlap)/ 3149 (staggered) configurations. For a set 
of 1102 configurations, we also measured the 12 lowest 
eigenmodes of both staggered and overlap operator. In 
all cases the quark mass was set to zero. 



The eigenvalues are ordered according to their imag- 
inary parts and the corresponding eigenvalue densities^ 
are plotted in Fig. [T] They display a gap-like behavior 
with the eigenvalue density starting to differ from zero 
considerably at aX ~ 0.15 and aX ~ 0.5, respectively. 
In addition, the overlap operator possesses exact zero 
modes, which we use to determine the topological charge 
of the configuration. 

To describe the RMT vs. Poissonian behavior of the 
Dirac spectra we will scan windows in the range of avail- 
able eigenvalues and measure the level spacing distribu- 
tions P{s) on unfolded eigenvalues [3], a typical quantity 
describing the eigenvalue statistics. The Gaussian RMT 
ensembles provide predictions for it, which also apply to 
systems with chiral symmetry, depending only on the uni- 
versality class. For gauge group SU (2) the latter are the 
Gaussian Orthogonal Ensemble (GOE) for the overlap 
operator (like for the continuum Dirac operator) and the 
Gaussian Symplectic Ensemble (GSE) for the staggered 
operator. The main difference of those P{s) formulae lies 
in the different repulsion strength of nearby eigenvalues, 
which results in a linear and quartic behavior of P{s) 
near s = 0, respectively. Independent eigenvalues, on the 
other hand, lack such a repulsion and the unfolded level 
spacing is a Poissonian distribution, i.e. P{s) = exp(— s). 

Fig. [2] shows the level spacing distribution of overlap 
eigenvalues (see [TU] for staggered spectra). It clearly re- 
veals that the level spacing agrees with the associated 
RMT predictions in the bulk and moves towards Poisso- 
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FIG. 1: Logarithm of the spectral density along the imaginary 
axis for the staggered (red, at smaller A) and overlap (blue) 
Dirac operator from the lowest 256 eigenvalues. 



^ Due to chirality, the nonzero eigenvalues come in pairs of opposite 
imaginary part and we restrict ourselves to the half with posi- 
tive imaginary part. In other words, all plots can be extended 
symmetrically around A = 0. 
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nian when the spectral window is shifted towards lower 
eigenvalues. To observe such a level spacing at least a 
few independently (Poissonian) distributed eigenvalues 
are needed on each configuration. Since independent 
modes occur only at the very low end of the spectrum 
where the spectral density is low, large enough volumes 
are required for that. In [5] also the independence of 
these data of the lattice spacing has been demonstrated 
for the staggered case. 

The properties of the independent and localized modes 
at the lower ends of the spectra are the main subject of 
the rest of the paper. We therefore check first, whether 
the modes of the overlap and the staggered operator see 
similar physical effects, meaning that their profiles are 
correlated and localized to similar locations. 

First of all we remark that for every overlap zero mode 
we find a staggered eigenvalue with unusually small value. 
In the topological sector Q = 0, the average smallest 
eigenvalue is 0.175(1), whereas it is 0.109(3) for |Q| = 1 
(and 0.098(8)/0. 141(6) for \Q\ = 2, where we have two 
small eigenvalues). This can also be seen in Fig. [7] bot- 
tom. 

0.07 r 
0.06 - 

0.05 - . -.^ ^ . .• 
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FIG. 3: Scatter plot of staggered (horizontal) vs. overlap (ver- 
tical) amplitudes for the lowest modes, I'i/jf (a;)| vs. |V'i^(a;)|, 
on a Q = configuration over the whole lattice. Data points 
with small amplitudes - as indicated by the circle - have been 
excluded to avoid overcrowding the plot. 



Next, the scatter plot of Fig. [3] gives a strong indica- 
tion for the correlation of the local mode amplitudes in a 
typical example configuration, which is further visualized 
by the two-dimensional profiles in Fig. [4j 

In order to quantify the similarity and localization of 
two modes we propose the following "interlocalization" 

i:^Nj2\rr:{x)\'\r:{x)\\ (i) 

X 

where N is the number of lattice sites, and 

is the absolute value (L^-norm) of the mth over- 




FIG. 4: Profile of the overlap (top) and staggered (bottom) 
lowest mode of the configuration of Fig. [S] in a lattice plane 
where the overlap mode takes on its maximum. (The absolute 
maximum of the staggered mode is separated from the overlap 
one by V2 lattice spacings in the remaining directions.) 
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FIG. 5: The interlocalization /, Eq. |T|, for matched modes 
(see text) as a function of the averaged overlap eigenvalue 
(aXm) (in lattice units), on a logarithmic scale and ensemble 
averaged. The maximal possible value for / on our lattice 
is 24^ • 4 ~ 5.5 ■ 10^, whereas delocalized modes yield / ~ 
1, indicated by the dashed gray line. Horizontal error bars 
visualize the spreads of the eigenvalues. 



lap/staggered eigenmode summed up over gauge - and 
in the case of overlap also spinor - indices at lattice site 
X. This is a positive quantity that receives large con- 
tributions when both modes are considerably large at 
some locations. Moreover, for two identical modes it be- 
comes their inverse participation ratio (IPR). The latter 
is a well-known measure for the localization, taking on a 
value of N for modes localized on a single point (on the 
lattice) and 1 for constant modes. In fact, we find J ~ 1 
also for two normalized modes with independent Gaus- 
sian distributed amplitudes at each site. Only modes 
that are similar and localized generate large values of /. 

We utilize / for matching the overlap and staggered 
modes^. We start by taking the lowest overlap mode 
and pair it with the staggered mode that has the largest 
interlocalization with it. Going up in the overlap spec- 
trum we continue this matching procedure in the same 
way, but use only those staggered modes that have not 
yet been paired up with a lower overlap mode. In Fig. [s] 
we plot the interlocalization values for the lowest modes 
matched in this way as a function of the corresponding 
overlap eigenvalue A™. From this plot it is clear that the 
zero modes (of the overlap operator, near-zero modes of 
the staggered operator) are matching close to perfectly. 
The value of 300 is actually in the same order of magni- 
tude as the IPR of the individual overlap zero modes^. 
Then / drops quickly and after a few modes it reaches 
the reference value 1 discussed above. 



^ In 1281 the positions of the highest peaks were used to reveal 
similarities between overlap and staggered modes. 

* The associated low-lying staggered modes have a slightly larger 
IPR, around 400, presumably because in contrast to the overlap 
operator the staggered operator is ultralocal. 



III. POLYAKOV LOOPS AS DEFECTS/TRAPS 

As an appetizer of our main finding we show in Fig. [6] 
the Polyakov loops (of one example configuration) in the 
lattice plane, where the lowest overlap and staggered 
eigenmodes take on their maximum, as shown in Fig. [4] 
The Polyakov loop is dominated by UV fluctuations at 
the scale of the lattice spacing as almost every lattice ob- 
servable. Therefore it is virtually impossible to see any 
structures in it. 

We applied 6 sweeps of APE smearing [291 EO] with 
a = 0.55 to the configuration, which leads to a smoother 
Polyakov loop landscape. Indeed, an island of Polyakov 
loop with opposite sign emerges at the location of the 
maximum of the lowest Dirac mode on the original un- 
smeared configuration. A similar profile becomes visible 
after simply averaging the (traced) Polyakov loops with 
their neighbors. 

Let us stress, that the lower panel of Fig. |6] is the 
only occasion that we present a smeared result. We now 
return to correlation functions of unsmeared Polyakov 
loops for the rest of this paper. 

To check the correlation of Polyakov loop islands 
and low Dirac modes in a quantitative way, we define 
"Polyakov loops as averaged by a particular mode", i.e. 
Polyakov loops weighted with the density of a normalized 




FIG. 6: Profile of the unsmeared (top) and smeared (bottom) 
Polyakov loop for the same configuration and in the same 
plane as in Fig. [4] By the naked eye nothing seems particular 
in this plane for the unsmeared case, whereas the smeared 
Polyakov loop actually takes its minimum (-0.68, compared 
to an average smeared Polyakov loop of 0.80) at the hotspot 
visible in the fermion modes in Fig. |4] (A correlation of the 
unsmeared Polyakov loop to fermion modes is exposed by 
virtue of statistical measurements, see text.) 
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Dirac mode, cf. [TB] . 

Lm ■= ^ \'ll)m{x)\'^ L{x) . 



(2) 



This quantity is restricted to the interval [—1, 1], just hke 
the Polyakov loop L. It is clear that a wave-like mode 
1pm with approximately constant amplitude yields an Lm 
close to the average Polyakov loop L{x)/V, whereas 
a strongly localized mode picks the Polyakov loop in that 
region. If the latter happens to be an island of "wrong" 
Polyakov loops, L„i tends to zero or even becomes nega- 
tive. 

Fig. [7] shows the ratio of averaged over differ- 
ent configurations and the average Polyakov loop, both 
for the low-lying overlap and low-lying staggered modes. 
The Polyakov loops averaged by the low-lying modes are 
indeed smaller than on average. Higher modes, on the 
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FIG. 8: The "Polyakov loop distribution as seen by the lowest 
overlap mode" , pi (L) from Eq. ([sjl , in the Q = sector com- 
pared to the Polyakov loop distribution and the Haar measure 
(valid in the low temperature phase). 



contrary, tend to see the average Polyakov loop. 

The connection between local Polyakov loops and low- 
lying Dirac modes is confirmed from a slightly different 
perspective by the following "Polyakov loop distribution 
as seen by a mode" . For that we weight the probability 
of Polyakov loops L by the amplitudes of the low-lying 
modes at those positions x where L{x) = L, 



P, 



(3) 



in continuous notion with some smeared delta-function 
Se (in practice we rescale histograms). The quantity Lm 
is recovered by the L-expectation value with this proba- 
bility. 



I 



dLpyn{L)L, 



( / dLpm{L) = 1) . (4) 



Thus the probability Pm visualizes how the global quan- 
tity Lm is generated by modifying the distributions of 
the local Polyakov loops. 

In Fig. [8] we show this probability for the lowest modes 
of the overlap operator, pi{L), on Q = configurations. 
One can clearly see that low fermion modes enhance low 
Polyakov loops down to L ~ — 1. 
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FIG. 7: The ratio of "Polyakov loops as averaged by low- 
lying modes" L,n, Eq. ([2|, to the average Polyakov loop for 
overlap (top) and staggered modes (bottom) as a function 
of the corresponding averaged eigenvalue (aXm). Horizontal 
error bars visualize the spreads of the eigenvalues. 



IV. INTERPRETATIONS 

A. Effective Matsubara frequencies 

In a background of constant temporal and vanishing 
spatial gauge fields the influence of the Polyakov loop on 
the Dirac spectra is very clear. We first diagonalize the 
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Polyakov loop introducing its phase 

Nt 



Y[ Uoixo) = exp {itp (^ _^ 



Xo = l 



(5) 



L = cos tp f £ [0,n] . (6) 



We gauge it into the last time slice, where it effectively 
changes the temporal boundary condition. 

The lowest modes of the free Dirac operator in the 
presence of this Polyakov loop are constant in space and 
plane waves in time, exp(ipxoT'). The quantum numbers 
p are governed by a combination of the antiperiodic tem- 
poral boundary condition for fermions and the Polyakov 
loop phase, namely p — t: ± tp + 27rZ, where the differ- 
ent signs emerge from the different color components, see 
Eq. ^. 

The eigenvalues of the free Dirac operator are these 
numbers multiplied by the temperature. The lowest ones. 



(7) 



we name effective Matsubara frequencies. These hold in 
the limit iVj — >■ c», whereas on the lattice one has 



1 



sm 



TT — (f 

Nt 



(8) 



At high temperatures the Polyakov loop (at fixed lat- 
tice spacing) becomes trivial, L — > 1, hence iy9 — > 
and the Matsubara frequency is Am — ^ sin (Tr/A^t). A 
more realistic estimate is obtained by using the average 
Polyakov loop at our temperature (L) — 0.37, from which 
we obtain the effective Matsubara frequency in lattice 
units (T = 1/A^ta) as 



Amo = si 



TT — arccos 0.37 



= 0.47, 



(9) 



which is the same order of magnitude as the lower end 
of the bulk of eigenvalues we measured (consistent with 
the findings of [H]). 

The main point of these considerations is that "wrong" 
Polyakov loops, L = —1 with ip = would lead to a 
vanishing effective Matsubara frequency. Am = 0, and 
thus to the lowest Dirac eigenvalues. 

Of course, in realistic configurations one has to take 
into account that the Polyakov loop varies in space and 
that nontrivial spatial links are present. Both will change 
the Dirac eigenvalues away from the free ones. Nonethe- 
less, the tendency that "wrong" Polyakov loops give rise 
to smaller eigenvalues persists and explains our finding 
about their pinning nature. 



Dirac modes in Yang-Mills theory (and QCD). The chi- 
ral condensate at zero temperature is thought of as due 
to instantons of realistic size and density and the first 
conjecture about the metal-insulator transition at finite 
temperature was based on instanton ensembles [3TJ [32] . 

The natural topological excitations at finite tempera- 
ture are magnetic monopoles. They appear as selfdual 
or antiselfdual solutions of the Yang-Mills equations of 
motion at finite temperature. As they are also elec- 
trically charged, they are called dyons. Dyons can be 
constituents of calorons (finite temperature instantons) 
[321 131] ! but may exist in isolation as well [35j . 

In gauge group SU{2), there are two dyons (and two 
antidyons) with opposite magnetic charge. One sort of 
these dyons is characterized by properties that exactly 
match our findings: the Polyakov loop at their core is 
— 1 and they possess zero modes with the physical an- 
tiperiodic boundary conditions [361137) (just like the con- 
stant configurations discussed above). The other sort has 
Polyakov loop and zero modes with periodic bound- 
ary conditions. 

In dilute ensembles of dyons most of the zero modes 
should remain low-lying modes. Then the first sort of 
dyons could explain the localized modes we analyzed, 
while the second sort could be responsible for low-lying 
periodic modes. At a positive average Polyakov loop one 
expects fluctuations to —1 much less frequent than those 
to -1-1, cf. Fig.[9j Hence the first sort of dyons could give 
independent low-lying modes, while the second sort could 
yield a condensate at periodic boundary conditions. This 
different appearance can be made quantitative by the dif- 
ferent fractions of unit topological charge the two dyon 
sorts have, which are such that the —1 dyons are indeed 
heavier as they have a larger (classical) action, see [10] 
for lattice evidence of this picture and [H] for a simple 
model. 

From the relation of magnetic monopoles to both low- 
lying modes and topological charge a crucial test of 
this picture is to compare these two quantities at given 
temperature and volume. Above the finite temperature 
transition fluctuations of the topological charge decrease 
sharply. This can be clearly seen by looking at how 
the topological susceptibility decreases at higher tem- 
peratures. As a result, at high temperature, topological 
objects presumably form a dilute gas of non-interacting 
objects. From the index of the overlap Dirac operator 
we have full information about the fluctuations of the to- 
tal topological charge. Assuming that in the dilute gas 
topological objects are uncorrelated, this can be used to 
compute the density of topological objects that in turn 
can be compared to the density of localized Poissonian 
Dirac eigenmodes. 

In Table [ij we show the probability of different topo- 



B. Topological objects 



Topological excitations of the gauge field and their zero 
modes are an attractive hypothesis to explain low-lying 



^ Monopoles also appear as defects |38j of the Polyakov gauge |39| , 
where they have L = —1 or L = +1 by definition. 
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2 




0.958(6) 


0.042(6) 


0.89(1) 


0.105(9) 


0.009(3) 



TABLE I: The probability of different charge sectors in the 
16'' X 4 (first two columns) and 24'^ x 4 ensemble (last three 
columns). The probabilities of charge sectors with the same 
magnitude but opposite sign have been added. 



logical sectors extracted from the index of the overlap 
Dirac operator. From the low occurrence of the topologi- 
cal charge sector ±1 and in particular of the sector ±2 it 
is indeed clear that topological objects form a very dilute 
gas. In these volumes it very rarely happens that there is 
more than one topological object on any given configura- 
tion. Since we need only an order of magnitude estimate 
of the density of topological objects we will ignore the 
probability of two or more objects occurring on any sin- 
gle configuration. In this approximation we can shortcut 
the calculation of the density of topological objects, and 
the average number of topological objects per configura- 
tion is just given by the probability of the \Q\ — 1 sector. 
On our volumes it is between 0.04 and 0.1 which is clearly 
far too small to account for the few localized Poissonian 
modes we found on average per configuration. 

This comparison rules out models based on uncor- 
related gluonic objects that carry both (0(1)) topo- 
logical charge and (antiperiodic) zero modes. How- 
ever, it does not exclude combinations of topological 
objects in which the topological charge cancels, like 
instanton-antiinstanton molecules originally suggested to 
be present with light dynamical quarks j42j or molecules 
of dyons called 'bions' carrying one near zero mode 
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FIG. 9: Independence of local Polyakov loops. Top: 
the auto-correlation A := {{L{x)L{y)} — {L{x)) {L(y))) / 
((L^(a;)) — {L{x))'^^ as a function of the lattice distance 
\x--y\/a (containing the same information as the free energy). 
Bottom: the level spacing distribution (of the Polyakov loop 
trace L) compared to the Poisson distribution. 



V. RANDOM MATRIX MODEL USING THE 
STAGGERED DIRAC OPERATOR 

Transitions between correlated and uncorrelated eigen- 
values like the one observed in the lattice data can be de- 
scribed by RMT models in various ways. One of the sim- 
plest possible ansatze is to start with a diagonal matrix, 
whose entries are uncorrelated random numbers, and add 
a matrix taken from one of the Gaussian ensembles. This 
model shows the desired transition in between parts of 
its spectrum with different eigenvalue density, however, 
the interpolating spacing distributions are different from 
the ones in the spectrum of the staggered operator [ilj . 
A possible explanation is the sparseness of this opera- 
tor and the spatial information that is contained in the 
next-neighbor interactions. In contrast, the full matrices 
from the Gaussian ensembles blindly connect all diagonal 
elements with equal strength. 



A. Motivation 



We construct a better suited random matrix model, 
based on sparse matrices, in the following. This model 
can be nicely motivated by our previous findings. Con- 
cerning the Polyakov loop, which will be connected to a 
random potential, three properties are relevant: (i) the 
distribution of local Polyakov loops extends to negative 
values and thus small effective Matsubara frequencies, cf. 
Fig. |8j (ii) Polyakov loops become independent quickly 
with their distance and therefore (iii) the level spacings 
of the Polyakov loop trace L - neglecting the spatial in- 
formation - are distributed according to a Poissonian dis- 
tribution. The last two properties are depicted in Fig. [9] 

Hence the Polyakov loops could provide the Poissonian 
ingredient for the statistics of the Dirac eigenvalues. To 
become more concrete, we remind the reader of the defi- 
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FIG. 10: Spacing distributions in several parts of our RMT spectrum plotted along with the GSE prediction and the Poissonian 
distribution. The parameters were fixed as discussed in the text. The data was obtained by an ensemble average over 2000 
random matrices. 



nition of the staggered Dirac operator, 

Dxx' = ^ X! [5x+i,,x-U^{x) - 5x-i,,x'Ul{x)\ , 



(10) 

with ri^{x) = (-1)2^-<M^" and U G SU{2). We split this 
operator in the temporal and spatial part 



(11) 



D^^ contains the hopping terms in the temporal direc- 
tion, whereas spatial hoppings are included in D^^ . 

The temporal part has a block-diagonal structure, with 
one block for each spatial site. We want to approximate 
each block and thus the whole Dirac operator by restrict- 
ing it to the subspace of the smallest eigenvalue quadru- 
plet of the temporal operator at each spatial lattice site. 
The quadruplet consists of two eigenvalue pairs with op- 
posite sign. The plus-minus degeneracy is necessary to 
conserve chiral symmetry, whereas the exact two-fold de- 
generacy^ has to be kept in order to have the right RMT 
universality class we find in the SU(2) staggered spectra. 

In this basis, the temporal part of the Dirac operator 



can be brought in the form 



with 



/- 


-Qx 





o\ 
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Ox) 



(12) 



1 . n - ips 
a \ Nt 



(13) 



which resemble the effective Matsubara frequencies on 
the lattice, as given in Eq. ([8|, however as a function 
of the local rather than averaged Polyakov loop phase 
(y9j. Therefore, £)'^^'("=o) is a diagonal matrix with very 
weakly coupled entries, as argued above. Physically 
speaking, the temporal part represents a (chiral) random 
potential. 

The spatial part becomes in the restricted basis 



„SP(n=0) 



U-x^x-^-i '^x.x~\-i 
'^x,x-\-i Ux^x~\-i . 



(14) 



In RMT language, this is Kramers degeneracy of the GSE. 



where the 2-by-2 matrices U and V can be shown to 
represent real quaternions. 

These have a concrete meaning: U connects eigenvalue 
pairs of equal sign, and is therefore responsible for the 



9 



250 r 



200- 



150 



where the diagonal entries 



0.06 



1.1r 



1 - 

0.9 
A 0.8 

V 

> 0.7- 
0.6- 
0.5 
0.4: 



0.08 



0.1 , 0.12 



0.14 



0.16 



= t{TT - 4>g) 



(17) 



are random numbers constructed with an overall scale t 
and a random angle (p^ £ [0, tt]. These quantities are 
equivalent to the effective Matsubara frequencies dg in 
the continuum, Eq. Q , the temperature T and the angle 
if^ of the local Polyakov loop, respectively. 

For we have taken the empirical distribution of the 
angle arccosI/(x) of local Polyakov loops, by converting 
the fine histogram in Fig. |8] accordingly. This yields an 
asymmetric distribution of (j)^ between and tt with a 
maximum below tt/2 (i.e. at positive Polyakov loop). Its 
most important feature, however, seems to be the tail 
towards the "trapping" (j)g ~ n (negative Polyakov loop 
locally), since we have observed that most of the features 
of the model persist when using the Haar measure sin^ (p 
as distribution [not shown]. 

For the spatial part, one first of all needs to fix the 
periodicity of the underlying space, i.e. an integer Ng 
such that X is identified with x + N^i, with unit vectors 
i in each of the spatial directions. M^^ , like ]jSP{ri=o) ^ 
has nonvanishing entries only at positions that connect 
next neighbors in that space. Its blocks 



0.06 0.07 0.08 0.09 0.1 , 0.11 0.12 0.13 0.14 
OX > 



, '^x,x-\-i '^x,x-\-i I 



(18) 



FIG. 11: Top: Inverse participation ratios of the 25 lowest- 
lying eigenmodes for the RMT model, plotted versus the av- 
erage eigenvalue of that mode. The data was obtained by an 
ensemble average over 2000 random matrices. Bottom: diag- 
onal entries "as avera ged by the lowest modes" of the RMT 
model, i9m from Eq. (191, divided by the average diagonal 
entry, to be compared with Fig. [7] Horizontal error bars vi- 
sualize the spreads of the eigenvalues. 



GSE-like level repulsion between nearest neighbors. V 
on the other hand generates the gap around zero in the 
spectrum as it connects eigenvalues of different signs. 



B. Explicit construction of the model 

We propose a random matrix model based on matrices 

M = M^^ + M^P, (15) 

that consist of two parts. M^^ has the same diagonal 
structure like £)^^("=°), i.e. 

o\ 

-ds 
i?s 
V ds) 



consist of random real quaternions u and w, that are RMT 
counterparts of U and V from Eq. (14). We have taken 



TE 



Mis 

XX 



(16) 



them to be Gaussian distributed around zero with their 
mean deviations au,<Jy as parameters of the model. 

By rescaling all the random matrices M it is clear, that 
only the ratios of the scales i, cr^ and (t„ are relevant 
parameters. To fix them we have measured the ratio of 
the average determinants'" of the empirical quaternionic 
hopping terms finding {det V) / {det U) « 1.6^ and took 
this ratio over for the ratio of a'^/a'^. The remaining ratio 
cr„/t = 0.2 {(Ty/t = 0.32) was put in by hand to obtain 
desired properties of the RMT model, namely a gap at 
zero and a transition between a Poissonian and a GSE 
spacing distribution. 

The eigenvalue density and spacing distribution of this 
model are plotted in Fig. 10 for a spatial extent Ng — 12 
and t = 1/4. We observe a similar transition like in the 
spectrum of the staggered Dirac operator. Another fea- 
ture that this model shares with the staggered operator 
is the increasing localization of eigenmodes as the corre- 
sponding eigenvalues decrease, quantified by the IPR's in 
Fig. [TT] top. 

In order to measure the correlation of the lowest eigen- 
modes of this RMT model to the diagonal entries, we 



^ For real quaternions, the determinant is just the sum over all 
squared components. 
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again define the latter "as averaged by a particular mode" 
^?n^ (19) 

X 

One can see in Fig. [TT] bottom, that the low modes are 
indeed localized to "islands" of low which are equiva- 
lent to low Polyakov loops. Hence this important effect 
is shared by our random matrix model, too. 

VI. CONCLUSIONS 

In the present paper we found a possible explanation 
for the emergence of localized Poissonian modes at the 
low end of the high temperature QCD Dirac spectrum. 
We showed that the localized modes are strongly corre- 
lated with large fluctuations of the Polyakov loop. This 
lowers the effective Matsubara frequencies for modes con- 
centrated there. We argued that the lowest part of the 
Dirac spectrum consists of this type of eigenmodes. We 
verified this picture for eigenmodes of both the staggered 
and the overlap Dirac operator. As a side result we 
also demonstrated that the spatial structure of the lowest 
overlap and staggered modes is highly correlated. This 
shows that different discretizations of the Dirac operator 
are sensitive to the same type of gauge field fluctuations. 

We also looked at the topological charge fluctuations 
as given by the index of the overlap. Assuming that at 
high temperature topological objects form a dilute gas 
and are uncorrelated, we could safely rule them out in 
creating the localized low modes. 

Finally, we proposed a dimensionally reduced random 
matrix model. It is based on sparse matrices encoding the 
three dimensional nature of the problem through nearest 



neighbor couplings from a lattice Dirac operator. To the 
best of our knowledge this is an example of a new kind 
of random matrix models for QCD, where so far only 
full matrices have been used. Beside chirality and the 
spectral gap our model reproduces the transition from 
localized Poissonian to delocalized random matrix type 
modes observed in the lattice Dirac spectrum as well as 
the correlation of the localized modes to "islands" of low 
on site "potential" . 

It is instructive to compare the Dirac operator to the 
Hamiltonian of Anderson type models, see [151 |35] for 
reviews. In the latter case usually diagonal (on site) dis- 
order is responsible for creating the transition to localized 
eigenmodes. In the case of the Dirac operator, the on site 
terms do not seem to be relevant. In fact, in the staggered 
operator they are exactly zero. However, in our dimen- 
sionally reduced three dimensional effective model the 
non-zero fluctuating on site terms are dynamically gen- 
erated by fluctuations of the local Matsubara frequency 
resulting from fluctuations of the Polyakov loop. In this 
way our dimensionally reduced effective random matrix 
model is analogous to the Anderson model. It would be 
interesting to study further how the presence and details 
of the transition depend on the parameters and matrix 
size defining the random matrix ensemble. 
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